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VIBRATION  FREQUENCIES  OF  A  CIRCULAR  CYLINDER. 
OF  FINITE  LENGTH  IN  AN  INVISCID  FLUID 


Introduction 


The  customary  procedure  for  taking  into  account  the  effect  of  the 
surrounding  fluid  on  the  vibration  of  a  free-free  beam  is  to  increase  the 
magnitude  of  the  mass  per  unit  length  of  the  beam  by  the  added  mass  at  each 
section,  considered  as  a  two-dimensional  form,  with  a  correction  factor  for 

the  three-dimensionaMty  of  the  flow.  In  the  present  work  there  are  compared, 

\j_-  Tag, 

\  for  a  particular,  case?*  the  natural  frequencies  of  vibration  obtained  by  thj-o 

V,  ^  are  mrpaverf 

^strip-theory  1  procedure  with  the  values  given  by  a  more  exact  formulation 
of  the  problem,  in  which  the  effect  of  the  presence  of  the  fluid  is  incorpo¬ 
rated  properly  into  the  vibration  equation^^T^  For  this  purpose  the  flex¬ 
ural  vibration  of  a  uniform  circular  cylinder  of  finite  length  was  selected 
for  study. 

General  Procedure 


The  displacement  of  the  beam  in  a  mode  of  frequency  cO  will  be 
assumed  to  be  of  the  form 


t)  -  'ir(^)  -<Usn  Odt  (1) 

in  which  z  denotes  distance  along  the  axis  of  the  cylinder,  which  extends 
from  z  »  -1  to  z  =;  1  .  The  amplitude  of  the  vibration,  v(z)  ,  is  celled 
the  displacement  function  of  the  beam.  By  differentiating  (1)  with  respect 
to  time  one  obtains  CGv(z)  as  the  corresponding  velocity  function. 

The  velocity  potential  for  the  motion  of  the  fluid,  ,  may  simi¬ 
larly  be  expressed  in  the  form 

(a) 

in  terms  of  cylindrical  coordinates  with  origin  at  the  center  of  the  beam. 
Both  (j)  and  <p  satisfy  Laplace “s  equation. 

*  Numbers  in  [  ]  indicate  references  at  end  of  this  report. 

/ 


•2- 


The  maximum  kinetic  energy  of  the  beam  is  given  "by 

Tb  =  X  aajO  v  dj  ( s) 

where  a  is  the  radius  of  the  cylinder  and  jog  its  density*  assumed  to  be 
a  constant.,  For  flexural  vibrations  the  maximum  potential  energy  V  is  [2] 

v=ldElM^d  (4) 

■  where  E  is  the  modulus  of  elasticity*  Furthermore  it  will  be  seen  that  the 
amplitude  of  the  kinetic  energy  of  the  fluid  due  to  the  vibration  of  the  beam 
can  also  be  expressed  in  terms  of  the  displacement  function  in  the  form 

Tf  =  1  f  “"j'jfM  ,3)  ni)  4  4  (s) 

where  p^,  is  the  density  of  the  fluid*  Thus  the  amplitude  of  the  total 
kinetic  energy  is  expressible  in  the  form 

T  =  Ts+ Tf  =  j[i r  aaj>  (s) 

where  c?  (  ^  ^  }  is  the  Dirac  delta  function 

f't^4)dr'1  <7) 

If  the  "added-mass"  function  A(  ^  *  <|  )  were  known*  one  could  obtain 
the  vibration  frequencies  by  the  Rayleigb-Ritz  method  from  the  expressions  for 

V  and  T  in  (4)  and  (6).  Because  the  solution  of  the  potential- flow  problem 
is  obtained  from  an  integral  equation  which  may  be  solved  by  the  use  of  quad¬ 
rature  formulas*  it  was  convenient  to  approximate  the  integrals  for  T  and 

V  also  by  means  of  quadrature  formulas*  The  potential  end  kinetic  energies 
are  thus  expressed  as  quadratic  forms  in  a  finite  numbe^  of  terms  from  which 
the  natural  frequencies  of  the  body-fluid  system  can  be  obtained  as  the  eigen¬ 
values  of  the  potential- energy  matrix  with  respect  to  the  total-energy  matrix* 
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these  matrices  being  composed  of  the  coefficients  of  the  respective  quadratic 
forms  [1] . 

Quadratic  Forms  for  V  and  Ts 

An  unsuccessful  initial  attempt  to  solve  the  problem  will  first  be 

briefly  described .  A  quadrature  formula  was  used  to  replace  the  integral  for 

T  in  (l)o  Next  the  interpolation  polynomial  corresponding  to  the  quadrature 
s 

formula  was  differentiated  twice*  The  resulting  polynomial  was  squared  and 
a  quadrature  formula  was  then  used  to  approximate  the  second  integral*  This 
attempt  was  not  successful,  because  the  resulting  quadratic  form  for  the  poten¬ 
tial  energy  derived  from  the  second  integral  was  not  positive  definite.  It  is 
believed  that  the  difficulty  was  caused  by  the  numerical  approximation  of 
v"(j|-)  ,  the  second  derivative  of  v(^)  . 

To  avoid  the  numerical  differentiation,  it  was  decided  to  find  a 
suitable  interpolation  polynomial  for  v"  ( g)  and  then  by  integration  to  obtain 
an  expression  approximating  v(3-)  .  This  is  the  procedure  which  will  now  be 
described  in  detail. 

The  conditions  to  be  met  by  the  interpolation  polynomial,  p(g.)  , 
are  that  it  he  equal  to  v"  ( at  certain  prescribed  points,  ^  ,  and  also 
that  it  satisfy  the  same  end  conditions  as  the  function  v(^)  For  a  free- 
free  beam  -  that  is,  one  which  is  not  restrained  at  either  end  -  these  condi¬ 
tions  are 

=  0 

v"'gHr-  ti  =  0 

Therefore,  the  conditions  prescribed  for 

M  P<£)  =  ij  = 

(d)  Pat)  =  0 

(=)  P(±£)  •=  o 

If 

ftj*)  =  ^ 


(8) 

(9) 


p( a)  are 


dr 


where  the  prime  indicates  differentiation  with  respect  to  in  the  right 
member  and  omission  of  the  factor  with  k  »  j  in  the  left  one,  then  the  ex¬ 
pression  VU 

r  A. 

k  ‘S-3i'S ,'QrM  J 

satisfies  condition  (a)  because  at  z  a 

■n  <?■  <?i 

_  x. 

where  5,,  *  0  or  1  according  as  i  £  j  or  i  a  j  , 


For  an  interpolation  polynomial  which  does  not  include  the  end  points 

among  the  z.  ,  the  end  conditions  may  be  satisfied  by  including  a  factor 
2  2  /  ^^2  2 

(g  -  l)  /( ^  -  1)  which  reduces  to  unity  when  ^  a  ,  and  is  zero  when 

Z.  =  t  1  .  Also,  the  first  derivative  is  zero  when  5«  f  1  .  The  desired 
a  3  a 

interpolation  polynomial  will  then  be 


Pm=l  Siiilhl  d- 

*  V.  J 


do) 


The  next  step  is  to  determine  an  approximation,  q{ to  v(  2)  in 

terms  of  the  quantities  d.  »  This  may  be  done  by  taking _ 

J 

a  ^ 

V'Q)  =  =  \*  f  p(^)  dtt  clir  +  4-  Ca  (n) 

The  polynomial  p(u)  may  be  taken  in  the  form 


Pu) 


'Vt  ,  v: 

t <i-,) 


■ctj 


'a. 


j*1  <&-')*  if'(i,-sk)  « 


r  Eri  ^ 


p-i 


If  the  coefficients  of  the  polynomial  TT  are  known,  the 

coefficients  may  be  most  easily  found  by  synthetic  division  of  H  by  2^ 

Since  is  a  root  of  IT  ,  the  remainder  should  be  zero,  which  serves  as  a 


check.  For  instance,  if  the  ^  are  the  Gaussian  quadrature  points,  TT  will 
he  the  Legendre  polynomial  of  degree  n  with  a  leading  coefficient  of  unity. 

On  combining  ai  .1  constant  terms,  one  can  write  p(u)  as 


j=i-  p=i 


with 


Performing  the  integration  indicated  in  (11)  gives 


and 


with 


W=  Cd+Ci  d3 

I 


G:,  =  L  z.r„. 


p=i 


p  p  j 


and  with 


z, 


_  Vi 


P  +  5 


z  Vt 


Pt\3 


~r 


} 


ptl 


(12) 


ZL  ~  p+5  3  P+3  aP+l  v 

°i(^  "f-p,  ^  " (p+ixp+3)  ■*"  i&-,)h*c'}+c*  (i3) 


(14) 


f  Cp+A)(  p  +  5)  (p+2.)(pt3)  ?  CP  +0 

The  constants  and  Cg  are  functions  of  the  quantities  dj  . 

In  order  to  evaluate  them  it  will  be  assumed  that  the  bar  has  (a)  zero  linear 
momentum  and  (b)  zero  angular  momentum.  Condition  (a)  gives 


i  =  0 


Substituting  (13)  into  the  integral  and  evaluating  it  give 


(15) 
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where 


f  .  i 

•  j  ^  *»y '' 


i 


It  is  assumed  that  a  is  even. 

Using  condition  (b)  of  zero  angular  momentum  gives 


Ij  dZ  = 


0 


(16) 


and  in  a  similar  way  C.  ' is  found  to  be 


C  -  IL  &•  dj 

1  j=t  j  j 


where 


£*V  ((2m+7X2WfsXim+4)  (2m5X2m+3)(2mt2)  '  (zYr,fi)(zmi)Z1ri) 


(17) 


Thus  q(x^)  can  finally  be  expressed  as 

=  g  Q/^s 

with 

C‘i  = 

The  potential  energy  of  the  beam  (2)  can  now  be  expressed  as  a  quad¬ 
ratic  form  in  the  quantities  d^  .  Use  of  the  interpolation  polynomial  P(|-) 
to  approximate  v"(^)  gives  for  the  potential  energy 

v  =?  ^ir-  f  ^ 


gives 


Using  a  quadrature  formula  with  quadrature  points  ^  (  j  1,  2*  “qii) 

V  = 

since  P(^)  s=  .  The  quantities  are  the  weighting  factors,  of  the 
quadrature  formula.  In  matrix  form  V  can  be  written  as 

=■  dr A  d 


V 


(18) 
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with  d  a  column  matrix  composed  of  the  d,.  and  A  a  diagonal,  n  x  n  . 
matrix  composed  of  the  weights  of  the  quadrature  formula.  The  symbol  T  de¬ 
notes  the  transpose  of  the  matrix. 

The  kinetic  energy  of  the  beam  is  obtained  from  (3)  by  substituting 
q_(  for  v( -p  and  applying  a  quadrature  formula.  The  result  is 


I 


2  2  Yl/ 

irR  w  a 


E  |?L  r  Q:  d:  E  Clk  dk 

to  1  jir  tJ  J  to  Lk  K 


r»h 


with 


'YV 


Bik  -  r  R;  c,  Cu. 

This  expression  can  also  be  written  in  matrix  form  as 


■T,  =>  .  d1  g^d  ' 


T 


(19) 


where  ' 

Bs  =(Bjb)  . 

For  the  kinetic  energy  of  'the  fluid  we  have  [3] 

r  as 

'f  ■  Z  4“  art 

where  the  integration  extends  over  the  surface  of  the  cylinder.  In  terms  of 
cylindrical  coordinates  this  assumes  the  form 


J  ?Tt 


(20) 


The  potential  (j)  is  assumed  to  be  due  to  an  axial  distribution  of  doublets 
^(,(^-)  oriented  in  the  direction  9  aft/ 2.  Thus,  the  velocity  potential  is 


o 


where 

and 


R!  =  *a  +  /  +Q^)2  =  (HA* 


%  a,  A  COS  9 


=  Jty  d/YbB 


■8' 


At  r  *  a  we  obtain  then 

2  i  2,  2 

in  which  Rq  “  (  ^  “  S  )  +  a  ,  Also  the  boundary  condition  at  the  surface  of 

the  cylinder  yields 


-  03  & 

The  kinetic  energy  of  the  fluid  is  then 


The  kinetic  energy  can  hence  be  written  as 


m 


at  I 


hJzCL 


(22) 


Tr  =  M}b)  di)d$  (23) 


where 

Wl(-^  =  Ltj-^r+eej*1 

Since  m ( ^>  ^  )  is  small  except  in  the  neighborhood  of  ^  e>  ^  ,  where 
it  rises  to  a  sharp  peak,  it  is  desirable  to  reduce  the  rapid  change  in  the 
integrand  in  that  neighborhood  so  that  representation  by  a  quadrature  formula 
is  more  accurate.  This  may  be  done  by  writing  the  expression  as  [4] 


T 


=  -  I  f  ( -M;p)  ^  ^ .. 


F 


Substituting  quadrature  formulas  for  the  integrals  with  the  exception  of 


gives 


T. 


f  " 
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with. 


V: 


i  =  ^P  ,  Ai'=/^)  ,  Mij  =  Mfc'Si) 


and 


£i~| 

a^^-i)2ta1V'a- 


The  kinetic  energy  can  he  written  as 


*w  EV4 

=  J^£*R.pH..R.  +^(m£-n())r.i]  '  ’• '  '  <«> 


where 


In  matrix  form 


u;  =  L 

ft 


Tf  =  -  ■Tft.w.o?.  VT  W  /f; 


with 


VK-  > 

£  =( 
r  V/^ 


W  =  (Wij)  S  tjrt/ . . 

On  the  boundary  of  the  solid,  r  ®  a  ,  we  have  from  (22) 


=  i£. 

<?  d-1. 


A=a 


Hence,  substituting  for  (p  from  (21),  we  obtain 

w  ^  ^  =  [f^  [JrT-  ^  ) 


1=0 


and  thus 


(25) 


co  wij)  =-|W)(i  -  &)dt 


,  K 


(26) 


with 


vi  =  ,  Ai'=/^)  ,  Mij=  "%.$)) 


and 


ii~\ 

o^i-iy^osy^ 


The  kinetic  energy  can  he  written  as 


^Tn&+-K  Ikp)) 

=  £  * R;  Ein-i  Ry  +  ^  -n(  ))^  ]  («) 


where 


In  matrix  form 


u(  =  f 

a 


Tf  =  -  irrW  /± 


with 


v; 

TT=  ( TA 
\KI 


t±=hz 


W  =  (Wi;)  =  '{^  Mij  R;  +  R$j (m*  }) ;  r  I; ,.-...  Vb 

On  the  boundary  of  the  solid,  r  ®  a  ,  we  have  from  (22) 


A.  -  Gl* 


Hence,  substituting  for  (p  from  (21),  we  obtain 


't.sCb 


and  thus 


( 25) 


-o  «: 


(26) 
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This  is  the  integral  equation  for  the  unknown  distribution  of  doublets, 

,  t)  .  A  solution  to  it  will  now  be  approximated  by  use  of  a  quadra¬ 
ture  formula.. 

If  a  modification  in  the  integrand  is  introduced  as  before,  the 

result  is 


CO 


=  -{(j>ft) *-/<$)  k  +/Uj)jV^)^ 


with 


3  CL 


~  pWF  ppipT^ 

Substitution  of  a  quadrature  formula  for  the  first  integral  gives 
OJVl  =  + 

u  J  J  ^ 

Kij  =  te(SiAi) 

kL=Mji)  =  i '.Ml;  M  d-t; 


in  which 
and 


-3; 


I  - 


+ 


i  + 


HFfeti)i+a5^  +  U 


L±jc  I 
K+O'+aFf3? 


(27) 


The  relation  between  the  velocity  of  the  boundary  and  the  doublet 
distribution  can  be  written  more  conveniently  as 


where 


=  -£  Htj/J-i 

(H«)  =  (RiK;i  +  4i(te,'-8<) 


s>  =  £ 


and 
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and 


co  x r  =  ~  V\j^ 


(28) 


Thus. 


A 


»  -tdH  ' V 


(29) 


When  (29)  for  jl  is  substituted  into  (25)  for  the  kinetic  energy, 
✓ 

the  result  is 

*i  *>» 


But  from  (17)  we  have 


so  that  we  obtain 


Tf  (j’wH'r 

'  2  —  —  —  — 


djdj 


:  r,  = 

-  dT  e,  d 

a  * 


with 


Bp  *  C'WH"  C 


(30) 


Thus,  .in  summary,  for  the  potential  energy  we  have 

V  -  dLAd 

8 

and  for  the  kinetic  energy,  combining  T„  and  T  ,  we  have 

X  s 

T  =  Jr  a 

where 

B  =  Bs+  -S-Bf 

To  complete  the  analysis  it  is  necessary  to  compute  the  frequencies 
so  that  they  can  he .  compared  with  those  obtained  by  the  two-dimensional 


(18) 


(31) 
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assumption.  From  a  physical  view,  what  has  been  done  is  to  replace  the  eontin» 
uous  cylindrical  "beam,  a  system  with  an  infinite  number  of  degrees  of  freedom, 
by  a  system  with  n  degrees  of  freedom,  n  being  the  number  of  ordinates 
used  for  the  quadrature  formula.  The  partial  differential  equation  for  the 
beam  vibration  is  then  replaced  by  n  ordinary  differential  equations.  By 
treating  the  quantities  d^  as  the  amplitudes  of  a  system  of  generalized  coor» 
dinates,  one  obtains  the  frequency  determinant  equation  for  the  system  [1] 


where 


A  ~XB 


^  0 


(32) 


X  fe*  Afkigl  (33) 

E  Q?  , 

An  alternative  method  which  could  be  used  to  obtain  (32)  is  the  Rayleigh-Ritz 
procedure .  The  expressions  for  the  potential  and  kinetic  energies  have  been 
expressed  as  quadratic  forms  in  terns  of  n  unknown  quantities,  the  d^  . 
Since  the  sum  V  +  T  is  constant  for  a  free  vibration  of  the  system,  we  have 


djAcl  -  XdTBcl  =  <xnst.  (34) 

In  this  quadratic  form,  the  coefficients  are  to  be  chosen  so  that  the 

frequency  is  a  minimum.  Partial  differentiation  of  (32)  with  respect  to  each 
of  the  ch  yields  dynamical  equations  from  which  (32)  is  again  obtained. 

A  simple  example  will  be  worked  to  illustrate  the  method.  If  n  &  2  , 
then  the  Gaussian  quadrature  points  are  (“0.57735027,  0.57735027).  The  calcu¬ 
lations  for  this  case  follow 


E  = 


57735027.  .57735027' 


& 


XL 


(■trTwSFfc) 

t \-i 


i 
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F 


'  H250000  I,t250000\ 

1,4485716  / 


zif  -  r  { 


it 


UtL 


z 


<p«0(p^)  (p^Kpta.) 

.14338272  -.02616,9139  \ 

.1^938272  iozbuq&q ) 


p(p«) 


&  =  U)(E) 


&  = 


A2i8D48Q0  ..  117073/2 


,11707312  .21404800 


Vz 


■{; 


^  ti  ?M+s)(2mt4)(2r/fr3)  (2W+3j(2rtH2)(innv()  (zrm-i)2H)(iw-)j 

-.15535714  j  «  ;>2 


7) 


ei 


•0)2 

I 


=  -3LILi 


1 


\ 


V 


e-,= 


Z.  e-  ~ 

*-i  ^ 


(;I268I086  — ,  1 268108k  ) 

-.0732(4827  .073214827  \ 
\  .073214927  073214827  j 


G  i  j  -  Z{  Sj  +  fj  +  Glj 

c  - 


.00352347  ,03442081  ■ 

.03442021  -ooqsza^?  ^ 
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As  a  check  on  the  calculations,  the  frequencies  of  vibration  trim 
be  computed  for  a  beam  vibrating  in  vacuo.  In  this  case 


6<:j  -  ZI  Cfc.;.  CK,> 

J  fe;-i 

R*  =  \  ,  ^(,2 

B  =  CTC 


f\3.\Q(S\  -6.65615  \ 

V-6.fc5fc95  13 .Ib4<?  J 


Since  A  =  I_  the  problem  is  most  easily  solved  for  the  inverse  of 
the  eigenvalues  which  are  related  to  the  frequencies.  That  is. 


where 


Thus, 


|B  -  Ki|  =0 


1 3.106.9  xIO'-K 
-6.6,5695  xlO"'1 


-6,65695*  I0'4 

i3.i0fc9*ID-4-K 


A  convenient  method  of  solving  this  determinant  is  first  to  add 
the  second  column  to  the  first,  and  then  to  subtract  the  second  rovr  from  the 
first.  This  gives 


and.  then 


0 

6.4500*  10'  — K 


-n.7(63#xlo“VK. 

13. VO 69  *  |D'4-t< 


K,  =  6.4j5_0Oj<1£l1 


1550.39 


Kz-  19.7638*  itT^ 


X  v=  505 ,98 
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Comparable  values  are  available  from  Timoshenko  [2],  where  a  solu¬ 
tion  for  the  one-dimensional  flexural  vibration  of  a  free-free  beam  is  given - 
The  fourth  root  of  A  here  is  comparable  to  the  value  given  by  Timoshenko, 
Table  I  presents  a  comparison  of  fX  with  the  square  of  Timoshenko's  value. 
This  value  is  proportional  to  oo,  the  frequency  as  given  by  (55) 

Table  I  -  Comparison  of  Eigenvalues  for  Free-Free  Beam 
Vibrating  in  Vacuo 

Timoshenko  Quadrature  Method  (n  @  2) 

1st  mode  22,37  22.49 

2nd  mode  61.67  39.36 

To  include  the  effect  of  the  fluid,  it  is  first  necessary  to  calm- 
late  the  matrix  Bf  .  For  this  purpose  the  radius  was  taken  to  be  a  a  0.10. 
The  ratio  of  fluid  density  to  solid  density  used  was  0.365,  which  is  repre¬ 
sentative  of  aluminum  in  water. 
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A  comparison  with  values  obtained  using  two-dimensional  strip  theory 
can  be  made  by  dividing  the  values  in  Table  I  by  /l  +  o^j og  .  This  comparison 
is  given  in  Table  II. 


Table  II  -  Comparison  of  Eigenvalues  for  Free-Free  Circular 
Aluminum  Beam  of  Length-Diameter  Ratio  10  Vibrat¬ 
ing  in  Water 


Strip  Theory  : 

More  Exact  Theory 

Timoshenko  Quadrature  Method  (n  a  2) 

Quadrature  Method  (n  a  2) 

1st  mode 

19.15  19.25 

IS.  38 

2nd  mode 

52.78  33.70 

33.78 

Considering  that  the  system  has  been  allowed  only  two  degrees  of 
freedom  (n  »  2)  ,  one  sees  that  the  agreement  between  the  eigenvalues  for  the 
first  mode  is  remarkably  good.  In  this  mode  the  frequency  given  by  the  more 
exact  theory  exceeds  that  by  strip. theory  by  0.67  percent.  The  corresponding 
three-dimensional  correction  factor  obtained  from  a  tvo»node  vibration  of  a 
spheroid  of  length- diameter  ratio  10  [5]  is  1.09,  an  increase  in  frequency  of 
9  percent.  Calculations  of  the  eigenvalues  with  larger  values  of  n  are  pres¬ 
ently  under  way,  and  it  will  not  be  clear  whether  the  very  small  correction  to 
the  frequencies  given  by  strip  theory  is  a  valid  one  until  this  work  is  com¬ 
pleted. 
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